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ABSTRACT 


Using  the  Stansell  (1970)  elevation  correction  for  navigational 
satellites  together  with  the  double  pass  method  of  Anderle  (1971)  and 
Greely  (1971)  but  for  Individual  navigational  satellites  on  two  consecutive 
passes,  it  is  shown  that  uncertainties  in  longitude  position  can  be 
reduced  to  the  order  of  10  to  18  m for  land  stations  with  just  a few 
observations,  and  less  than  4 m for  dock-side  ship  stations  taken  over  a 
period  of  a week.  It  is  also  shown  that  the  method  allows  identification 
of  satellites  not  giving  reliable  data  and  permits  conclusions  to  be  reached 
as  to  the  cause  of  an  apparent  error  dependance  on  surface  elevation  noted 
above  1000  m in  the  Andean  region  of  Colombia,  South  America,  reported  by 
Woollard  and  Thompson  (1974).  An  additional  advantage  of  the  double  pass 
method  brought  out  is  that  it  allows  geoidal  heights  to  be  determined  for 
remote  islands  with  the  same  degree  of  reliability  as  their  position  and 
comparable  to  the  best,  continental  determinations  of  geoidal  height  using 
standard  geodetic  measurements. 
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INTRODUCTION 


Hawaii  Institute  of  Geophysics  Report  HIG-74-1  by  Daugherty  (Part  I) 
and  Woollard  and  Thompson  (Part  II)  deals  with  the  accuracy  of  satellite 
determined  positions  at  various  port  sites  around  the  Pacific  and  in  the 
Andean  region  of  Colombia,  South  America,  and  their  difference  from  map 
coordinate  positions. 

As  Daugherty  in  Part  I (hereafter  referred  to  as  Daugherty)  brings  out 
the  fact  that  successive  passes  by  different  satellites  with  the  ship  at  a 
fixed  position  for  periods  of  a few  days  to  two  weeks  indicated  that  erratic 
values  were  obtained,  there  is  some  question  as  to  how  much  value  can  be 
attached  to  a single  pass  satellite-determined  position  when  a ship  is 
underway.  Woollard  and  Thompson  in  Part  II  (hereafter  referred  to  as 
Woollard  and  Thompson)  bring  out  a different  type  of  problem;  namely 
systematic  differences  in  longitude  for  East- and  West-passing  satellites 
that  appeared  to  be  elevation-dependent  at  elevations  above  1000  m.  These 
data  have  been  reexamined  as  to  the  causes  of  the  discrepancies  noted. 

By  way  of  review,  Daugherty  used  a simple  arithmetic  averaging  of  the 
data  obtained  on  successive  passes  of  each  satellite  recorded  by  using  a 
MX/702/hp  satellite  navigation  system  on  University  of  Hawaii  ships  in 
1970,  1971  and  1972  at  eight  port  sites  distributed  throughout  the  central 
and  south  Pacific.  Daugherty  also  averaged  the  results  obtained  for  each 
individual  satellite  and  computed  the  standard  deviation  of  each  pass  from 
the  mean  as  well  as  that  of  the  mean  values. 

A summary  of  Daugherty's  results  with  the  number  of  observations  at 
each  site  which  varied  from  6 to  153  are  given  in  Table  1.  It  will  be  noted 
that  two  series  of  observations  were  made  at  Suva  with  significant  differences 
in  the  apparent  reliability  of  the  longitude  position  as  a result  of  an 
apparent  error  in  initialized  antenna  height  for  the  first  series.  In  the 
present  study  the  effect  of  the  error  (20  m)  is  evaluated  and  the  data  for 
Suva  used  to  demonstrate  the  effectiveness  of  the  double  pass  method  for 
determining  position  and  antenna  height  (geoidal  height)  at  sea  level  sites. 

Because  of  time  restrictions,  the  land  data  reported  by  Woollard  and 
Thompson  were  much  more  limited  in  terms  of  the  number  of  satellite  passes 
recorded  at  each  site.  The  apparent  ground  elevation  dependance  of  position, 
particularly  longitude  reported  for  East- and  West-passing  satellites , is 
brought  out  in  Table  2 (from  Woollard  and  Thompson,  their  Table  2) . 

Although  Woollard  and  Thompson  showed  that  if  an  empirical  correction 
for  the  tangent  of  the  satellite  elevation  angle  were  applied,  reasonable 
agreement  in  position  could  be  obtained  with  available  maps,  principally 
the  ONC  scale  1:1,000,000  series,  the  discrepancies  at  high  elevations  still 
exceeded  those  found  in  general  for  sea-level  sites. 

It  was  with  the  above  datum  sets  that  the  writer  started  the  investiga- 
tion here  reported.  As  part  of  the  investigation, the  reliability  of  maps 
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available  in  Colombia  was  also  investigated  since  different  map  series  for 
the  same  area  indicated  differences  in  coordinates  for  the  same  locations. 

The  principal  differences  in  the  writer's  approach  from  that  used  by 
Daugherty  and  that  used  by  Woollard  and  Thompson  was  in  the  application  of: 
(1)  the  Stansell  (1970)  correction  curve  for  antenna  height,  which  while 
similar  to  the  tangent  correction  incorporated  by  Woollard  and  Thompson  for 
their  corrected  data,  also  appears  to  incorporate  tropospheric  refraction 
effects,  and  (2)  the  double  pass  method  of  Anderle  (1971)  and  Greely  (1971), 
published  as  U.  S.  Naval  Weapons  Laboratories  technical  reports  that  came  to 
the  writer's  attention  in  an  article  by  Schultz  et  al.  (1975)  after  most  of 
the  work  for  this  report  was  completed. 

As  the  writer's  interest  in  this  problem  was  generated  through  his 
participation  in  Project  Narino.,  a seismic  crustal  study  in  Colombia  which 
was  depending  on  the  satellite  defined  positions  reported  by  Woollard  and 
Thompson  for  the  location  of  the  land  shot  and  recording  sites,  most  of  the 
report  is  devoted  to  this  series  of  measurements  and  the  reliability  of  the 
maps  used  in  connection  with  Project  Narino. 

Because  of  the  large  sample  of  data  available  for  Suva  and  the  fact 
that  two  different  series  of  observations  had  been  made  with  different 
initialized  antenna  heights,  these  data  were  also  studied  to  see  whether  the 
same  approach  as  applied  to  the  Colombia  data  would  significantly  reduce 
the  uncertainty  in  position  found  by  Daugherty.  As  will  be  shown,  this  was 
indeed  the  case,  and  from  the  study  it  became  apparent  that  the  satellite 
data  offers  a direct  method  for  determining  geoidal  height  having  the  same 
degree  of  reliability  as  that  of  the  derived  longitude  position. 

In  presenting  this  report  the  internal  agreement  between  the  two  sets 
of  maps  used  on  Project  Narino  will  be  discussed  first.  These  were  the 
Instituto  Geografico  "Agustin  Codazzi"  (IGAC)  maps  on  a scale  of  1:25,000 
and  maps  prepared  by  the  Texas  Oil  Company  (Texaco)  on  a scale  of  100,000. 
Both  were  based  on  the  Bogota  datum.  The  reading  accuracy  for  comparisons 
between  the  two  sets  of  maps  was  about  0.015'.  As  will  be  shown,  discrepan- 
cies between  the  two  sets  of  maps  appear  to  lie  within  the  accuracy  with 
which  positions  could  be  read  from  the  Texaco  maps.  On  the  whole,  agreement 
was  good  except  in  Buena  Ventura -Val le  area,where  a consistent  discrepancy 
of  about  0.6'  in  latitude  was  found  to  characterize  part  of  the  Texaco  map. 

The  second  section  is  a discussion  of  the  double  pass  method  for 
determining  longitude  and  height  using  the  Stansell  (1970)  elevation 
error  curve  for  satellite  elevation.  The  third  section  presents  the  results 
for  the  redetermination  of  positions  in  Colombia  using  the  Stansell  (1970) 
satellite  elevation  correction  and  the  double  pass  method.  In  the  fourth 
section  the  same  methodology  is  applied  to  the  large  body  of  data  available 
for  Suva.  As  will  be  seen, this  last  indicates  that  positioning  in  longitude 
and  height  is  better  than  8 m for  a daily  mean,  better  than  4 m for  a 
weekly  mean,  and  about  2 m if  a two- week  mean  is  used,  without  earth-  or 
ocean-tidal  corrections. 
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IGAC  AND  TEXACO  MAP  COMPARISON 


Both  the  IGAC  map  and  Texaco  map  coordinate  systems  are  based  on  a 
gauss  conformal  projection  with  the  origin  at  the  National  Observatory  in 
Bogota,  Colombia,  at  4°35'56"570N  and  47°04'51"300W  and  X = 1000N  and 
Y = 1000E.  Most  of  the  available  IGAC  maps  were  on  a 1:25,000  scale  and 
the  Texaco  maps  were  on  a 1:100,000  scale  (with  a 3°  E-W  offset).  A 
computer  program  used  here  to  transform  the  X,  Y coordinates  to  latitude 
and  longitude  is  based  on  J.  H.  Addison's  program  which  was  kindly  supplied 
by  J.  F.  Gettrust,  now  at  the  Australian  National  University  in  Canberra, 
formerly  at  University  of  Wisconsin,  and  adapted  for  the  University  of 
Hawaii  computer  by  A.  Kam.  The  constants  used  in  the  program  are: 

A = 6378388  (radius)  and  ES  = 0.00672267  (eccentricity^). 


Coordinate  transformation 


A check  was  provided  by  reading  the  x and  y coordinates  of  the  map 
corners  of  some  Texaco  maps  (where  \,  tp  were  also  given)  and  computing 
X,  cp  from  x,  y.  These  calculations  are  presented  in  Table  3.  The  C: 
line  gives  the  map  designation  (such  as  L-l)  and  the  geographical 
coordinates  in  degrees -minutes  as  indicated  on  the  map;  the  following  lines 
gives  the  x and  y coordinates  measured,  and  the  calculated  geographical 
coordinates  in  degrees -minutes , N latitude  and  W longitude.  The  standard 
deviation  from  the  indicated  coordinates  is  01013  in  the  N-S  and  010076  in 
the  E-W  direction;  these  values  include  the  measuring  error  and  also  the 
drafting  inaccuracies  (reduced  to  map:  0.24  mm  NS  and  0.14  mm  EW) . 

We  can  therefore  conclude  that  the  computer  program  does  not  contain 
any  gross  errors. 


Inter-map  comparison 

To  test  the  agreement  between  the  IGAC  and  Texaco  maps,  several  areas 
were  selected  for  comparison.  The  selection  was  based  on  the  proximity  of 
the  areas  to  the  seismic  refraction  lines,  satellite  navigation  fix  points, 
and  availability  of  maps.  In  each  area  a number  of  easily  identifiable 
features  were  used  for  comparison  points,  such  as  river-bends,  or  road-river 
or  railroad  intersections.  Comparisons  for  the  areas  are  given  in  Table  4. 
Each  point  comparison  consists  of  four  printout  lines: 

Line  1 : gives  the  area  (at  only  the  first  mention) 

BUENA  VENTURA -VALLE,  IGAC  map  designation  260-III-C, 
and  Texaco  map  designation  J-4; 
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Line  2:  IGAC  map  input  x,  y;  output  degrees,  minutes  N,  W; 

Line  3:  Texaco  input  x,  y;  output  degrees,  minutes  N,  W;  and 

Line  4:  difference  in  output 

degrees,  minutes  N,  W. 

The  difference  in  line  4 is  positive  for  the  IGAC  position  being  to 
the  north  and  to  the  west  of  the  Texaco  position. 

In  general,  the  use  of  small  river-bends  as  comparison  points  gave  the 
best  results.  On  the  Texaco  map,  roads  at  times  appeared  to  have  been 
smoothed  out;  or  just  one  side  of  the  road  located  in  its  correct  position 
with  the  other  drawn  parallel  to  it.  The  comparisons  were  relatively  poor 
in  the  southern  area  of  the  Buena  Ventura-Valle  area  since  the  available 
maps  contained  few  good  comparison  points.  However,  serious  discrepancies 
do  exist.  Some  part  of  the  NW  corner  of  the  J-4  Texaco  map  must  have  been 
shifted,  as  can  be  seen  by  the  016  southward  shift  for  the  north  part  and  a 
0.74  westward  shift  for  the  south  part  with  respect  to  the  IGAC  map.  The 
shift  is  clearly  marked  on  the  J-4  map,  where  a southbound  road  ends,  with  a "?" 
sign,  some  1 km  NE  of  the  bend  in  the  Buena  Ventura  road  (heading  EW) 
instead  of  joining  it  as  it  does  on  IGAC  map  260-III-C.  (Coordinates  of 
the  Buena  Ventura  road-bend  x = 921.78,  y = 678.68;  that  is  the  fifth  point 
in  the  Buena  Ventura  area  in  Table  4.) 

Table  5 summarizes  the  results  of  Table  4 by  area.  It  indicates  the 
N-S  and  E-W  offsets  for  each  area  and  the  standard  deviation  of  a single 
observation.  Most  of  the  offsets  noted  are  close  to  one's  ability  to 
measure  features  correctly  and  to  the  accuracy  of  the  drawing  of  the  Texaco 
map  as  reflected  by  the  standard  deviations  of  about  01015  (or  near  28  m) . 

Other  than  in  the  Buena  Ventura  area,  on  the  eastern  flank  of  the  Cordillera 
(Mocoa)  some  0!06  N-S  and  E-W  offsets  appear  that  are  larger  than  the 
standard  deviation  of  a single  observation  based  on  four  points. 

In  conclusion,  there  are  serious  offsets  (probably  internal  to  the 
Texaco  map)  in  the  Buena  Ventura  area  and  a small  one  west  of  Mocoa. 

Except  for  those  in  the  Buena  Ventura  area,  none  of  the  other  small 
discrepancies  should  affect  the  positions  of  seismic  refraction  stations 
with  respect  to  the  geographical  coordinates  from  either  set  of  map 
coordinates . 
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DOUBLE  PASS  METHOD  APPLICATION  TO  LONGITUDE  AND  HEIGHT 
DETERMINATION  OF  THE  SATELLITE  NAVIGATOR 
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Elevation  angle  correction 

Large  differences  (500  m or  more)  appeared  when  the  satellite 
navigation  positions  adopted  by  Woollard  and  Thompson  were  compared  with 
the  known  map  positions,  especially  for  Tumaco  and  Laguna  la  Cocha.  The 
question  arose  as  to  whether  the  tangent  correction  applied  was  valid. 

To  test  this,  the  data  for  Laguna  la  Cocha  longitudes  were  plotted  as  a 
function  of  the  tangent  of  the  elevation  angles  (Fig.  1).  As  seen,  the 
lines  fitted  to  the  E-passing  satellites  and  the  W-passing  satellites  do 
not  intercept  the  zero  elevation  line  at  a common  longitude.  The 
resulting  gap  is  of  the  order  of  900  m.  The  computer  program  used  by 
Woollard  and  Thompson  (taking  the  least  squares  fit  of  all  points) 
obviously  provided  a position  somewhere  in  the  gap  as  did  their  graphical 
solution  using  uncorrected  data.  The  tangent  correction,  while  clearly  a 
step  in  the  right  direction  on  the  basis  of  the  model  used  by  Woollard 
and  Thompson  to  explain  the  cause  of  the  discrepancies  observed  in  longi- 
tude position  as  being  an  error  in  antenna  height  incorporated  in  the 
satellite  software  program,  was  not  the  complete  answer.  The  answer 
itself  became  apparent  from  a paper  by  Stansell  (1970)  provided  the  writer 
by  Dr.  John  Rose. 

In  this  paper,  Stansell  (1970),  shows  a graph  (his  Fig.  11)  wh.ch 
gives  the  fix  error  (in  longitude)  as  a function  of  satellite  elevacion 
angle  K that  results  from  the  error  in  the  antenna  height  estimate 
(initialized  antenna  height).  If  f(K)  is  the  fix  error  (in  n.  miles) 
resulting  from  a 1-m  height  error,  then  a dH-meter  height  error  will 
result  in  a dH  f(K)  n.  mile  error  in  longitude. 

If  a number  of  satellite  passes  are  recorded,  each  having  elevation 
angle  and  fix  L^,  and  data  points  for  satellite  pass  no.  ^ are  plotted 
in  a rectangular  coordinate  system  at  points  y^  = f(K^),  x^  = L^  with  no 
other  noise  present,  then  th^  points  y^,  x^  will  fall  on  a straight  line 
for  which  the  slope  is  = 3h  * T*le  ^nverse  sl°Pe  will  give  the  antenna 
height  error  and  the  value  for  L at  y = 0 = f(K=0)  will  be  the  correct 
longitude  position. 

Since  the  input  data  were  not  known  and  no  computer  programs  were 
available  for  calculating  Stansell 's  (1970)  graph,  values  were  carefully 
measured  from  the  Stansell  diagram  (given  for  changes  in  L for  100-foot 
errors  in  antenna  height)  and  recalculated  for  the  position  changes  that 
would  result  from  a 1-m  antenna  height  change.  These  values  for  10° 
increments  of  elevation  angle  were  plotted  on  a large  graph  and  a smooth 
curve  was  drawn,  from  which  the  final  values  were  read.  (The  approximate 
values  for  f(K)  and  K are  listed  in  Table  6 and  also  shown  in  Figure  9.) 
Position  changes  that  would  result  from  a similar  correction  curve  given 


in  the  Magnavox  706  manual  are  slightly  smaller  than  the  ones  derived  from 
the  Stansell  curve.  Therefore,  antenna  height  changes  based  on  the 
Magnavox  curve  would  be  slightly  larger  than  the  ones  used  here. 

In  Figure  2,  a plot  of  f(K^)  versus  (the  longitude  fixes)  for 
Laguna  la  Cocha  is  shown  and  the  inverse  of  the  slope  defined  indicates 
an  elevation  error  of  some  2000  m.  Correcting  the  longitudes  for  an 
error  of  dH  = 2110  m and  plotting  again  the  reduced  versus  f(K^)  values 
results  in  a remaining  spread  of  some  0.2  min  or  nearly  0.2  n.  miles 
(neglecting  cos  (latitude))  (Fig.  3). 

As  the  direction  of  the  satellite  navigator -computed  longitude  is  the 
opposite  from  the  direction  of  the  passing  satellite  if  the  initialized 
antenna  height  is  too  high  (that  is  toward  the  west  for  an  east-passing 
satellite),  or  toward  the  satellite  if  the  initialized  height  is  too  low, 
it  is  clear  that  the  satellite  navigator  program  (for  Laguna  la  Cocha  and 
other  high -elevation  stations)  used  a different  (too  elevated)  antenna 
height  than  the  one  punched  in  by  the  equipment  operator  (N.  J.  Thompson). 
This  will  be  discussed  in  more  detail  later. 


Double  pass  method 


As  can  be  seen  from  Figure  3 the  spread  in  the  longitude  positions 
can  be  reduced  further  (to  0.01  min  or  20  m)  if  the  data  are  restricted  to 
height  corrected  values  for  individual  satellites  for  which  there  are  two 
consecutive  passes.  In  the  plot  the  value  of  f (K)  for  the  elevation  angle 
is  the  y coordinate  and  the  corrected  longitude  the  x coordinate  with 
E-passing  satellites  plotted  positive  for  f(K)  = 0 and  W-passing 
satellites  plotted  negative.  The  intersection  of  the  resulting  lines  of 
position  with  f(K)  = 0 defining  the  more  precise  position.  When  64  double 
passes  were  analyzed  for  Suva,  the  standard  deviation  of  a single  such 
determination  was  found  to  be  about  20  m in  longitude  and  about  the  same  in 
height  (Tables  12  and  13).  Part  of  this  standard  deviation  seems  to  be  due 
to  differences  in  individual  satellites,  notably  no.  64  in  1973  in  Colombia 
(the  same  as  no.  65  in  1971  in  Suva)  which  shows  a lower  antenna  height 
(Figs.  3,  16,  and  18).  In  passing  it  should  be  noted  that  in  general,  the 
longitude  position  data  for  the  double  passes  in  Colombia  in  1973  seem  to 
indicate  less  spread  than  those  recorded  on  board  the  KANA  KE0KI  in  Suva 
in  1971.  The  lower  accuracy  for  Suva  is  probably  due  to  effects  of  the 
ship's  motion  with  wave  and  tidal  action  in  the  harbor  and  perhaps  to  better 
satellite  orbital  parameter  programming  in  1973. 

Note  that  no  corrections  have  been  applied  for  the  tropospheric 
refraction  index,  for  which  the  satellite  navigator  computer  assumes  an 
average  value. 


COLOMBIA  SATELLITE  NAVIGATOR  POSITIONS 


Observation  polnts--IGAC  map  positions 

Woollard  and  Thompson  presented  the  satellite  navigation  data 
together  with  site  descriptions  and  compared  them  to  positions  on  the 
1 : 1 ,000 ,000-scale  ONC  L26  aeronautical  chart.  This  chart,  the  only  one 
covering  all  sites  on  the  same  projection  and  datum,  was,  however, 
inadequate  for  precise  determinations  of  location  because  of  the  scale. 
Therefore,  Thompson,  who  made  the  field  observations  in  Colombia,  and  the 
writer  reviewed  site  positions  where  possible  on  the  IGAC  maps  with 
scales  of  1:25,000,  1:10,000,  and  in  one  case,  1:2,000.  Longitudes  and 
latitudes  were  then  computed  from  the  x,  y coordinates  using  the  program 
supplied  by  G.  F.  Gettrust.  The  determined  positions  and  elevations  are 
given  in  Table  8.  However,  even  better  accuracy  could  be  achieved  if  the 
sites  could  be  resurveyed,  especially  the  four  for  which  photographs  of 
the  antenna  positions  are  available.  For  convenience,  the  original 
satellite  positions  and  heights  from  Woollard  and  Thompson  are  presented  in 
Table  7.  As  it  appears  that  the  site  elevations  (surface  elevation  plus 
geoid  height)  were  not  used  consistently  in  their  report,  the  antenna 
height  and  geoid  height  as  initialized  in  the  field  were  taken  directly 
from  Thompson's  original  book  and  added  as  two  additional  columns  in  Table  7. 

The  estimated  accuracies  of  x,  y,  and  height  as  presented  in  Table  8 
(and  based  on  available  photographs,  scale  of  the  maps  used,  and  discussions 
with  Thompson)  with  respect  to  the  maps  are  as  follows: 


Site 

(N-S)  x (in  m) 

(E-W)  y (in  m) 

Height  (m) 

Buga 

± 10 

± 20 

± 10 

Popayan  I 

± 2 

± 3 

± 3 

Laguna  la  Cocha 

± 30? 

± 20? 

± 20? 

Tumaco 

± 15 

± 15 

+ 2 

However,  an  accuracy  of  the  order  of  20  to  50  cm  could  be  obtained  if  the 
sites  were  to  be  surveyed  and  based  on  the  available  pictures  of  the  antenna 
pos it  ions . 


Magnavox  706  height  initialization  error 

The  satellite  longitude  fix  versus  elevation  angle  function  f(K)  for 
the  Colombia  positions  are  presented  in  Figures  2 through  7.  As  Figures  2 
and  3 indicate,  a very  large  antenna  height  correction  at  Laguna  la  Cocha 
is  required  to  reduce  the  data  by  the  double  pass  method  to  the  final  mean 
longitude  of  77°08'8050  ± 0’0033  (crM) . Since  Cali  (Fig.  5)  and  Popayan  I 
(Fig.  6)  also  show  large  height  corrections,  an  operator  error  for  the 
height  Initialization  is  excluded. 


8. 


The  best  explanation  that  appears  reasonable  for  the  capricious 
height  initialization  of  the  Magnavox  internal  computer  will  be  discussed 
relative  to  Laguna  la  Cocha.  As  shown  by  Figures  2 and  3,  the  computer- 
initialized  height  is  too  high  by  2110  + (+7-27-153) /3  = 2052  m.  This 
means  that  the  computer  was  working  with  an  initial  height  somewhere  above 
2800  m plus  2052  m,  or  over  4800  m,  if  the  site  altitude  of  2808  m is 
approximately  correct.  This  could  happen  as  follows  (borne  out  also  by  the 
Cali  and  Popayan  data):  the  computer  dislikes  altitudes  over  999  (or 

1024?)  m.  When  an  elevation  over  that  height  is  initialized,  the  computer 
doubles  the  corresponding  decimal  (1000)  or  binary  digit  (1024).  The  data 
are  too  few  and  the  uncertainty  in  the  absolute  value  of  the  Stansell 
correction  too  large  to  determine  which,  but  it  probably  is  the  binary. 

The  result  then  is  the  following: 

Operator  H = 2808  = 2048  + 760 
Computer  H = 2 x 2048  + 760  = 4856 
Concluded  H = Computer  H - height  correction 
= 4856  - 2052  = 2804  ± 60  tn  (aM) . 

Until  the  uncertainty  in  computer  shift  of  24  m/1000  m at  elevations  higher 
than  1000  m is  resolved,  the  height  data  for  Laguna  la  Cocha  (48-m  shift) 
and  Popayan  and  Cali  (24-m  shift)  should  not  be  used  without 

putting  these  possible  errors  in  evidence.  The  Bogota  data  are  inconsistent, 
probably  because  of  interference  problems,  and  have  not  been  considered 
reliable  for  the  double  pass  method. 


Colombia  satellite  positions — comparison  with  IGAC  coordinates 

The  satellite-determined  coordinates  and  the  IGAC  map  coordinates  are 
compared  in  Table  9.  In  this  table  the  satellite  latitudes  are  those 
determined  by  Woollard  and  Thompson  using  intersecting  lines  of  position 
for  each  satellite,  except  for  Tumaco,  for  which  an  arithmetic  mean  value 
was  used.  The  satellite  longitudes  are  the  values  obtained  here  from  the 
double  pass  method  and  are  derived  from  the  f(K)  "0"  crossings  in  Figures  2 
through  7.  The  cfy  errors  reflect  the  internal  consistency  of  the  data.  The 
purpose  of  analyzing  the  Suva  shipboard-fixed  station  data  in  the  next 
section  of  this  report  is  to  substantiate  the  small  errors  found  by  applying 
the  double  pass  method  in  particular  to  the  Colombia  data.  All  the  Colombia 
satellite-determined  longitudes  are  probably  within  10  m,  and  similarly  the 
heights  for  Buga,  Tumaco,  and  Buena  Ventura  probably  within  better  than  20  m. 
The  height  error  evaluation  for  Cali,  Popayan  I,  and  Laguna  la  Cocha  will 
depend  on  the  Magnavox  706  programming  error  and,  because  of  the  large 
corrections  involved,  on  the  accuracy  of  reading  the  Stansell  elevation 
angle  correction.  The  height  correction  for  Tumaco  was  obtained  by  best -fit 
combining  satellites  no.  41,  53,  and  63  in  pairs  and  determining  the  slope 
of  the  corresponding  lines  (Fig.  7).  Satellite  no.  64  was  not  used,  since 
height  determination  from  no.  64  seems  distinctly  different,  as  can  be  seen 
for  Laguna  la  Cocha  (Fig.  3)  and  the  corresponding  determinations  in  Suva 


(satellite  no.  65)  (see  Table  13  and  Figs.  16  and  18).  However,  longitude 
positions  for  no.  64  do  not  deviate  from  those  obtained  by  the  double  pass 
method  from  the  other  satellites.  The  no.  64  (In  1973)  Colombia  Is  the 
same  satellite  designated  no.  65  (In  1971)  In  Suva. 

So  far  no  IGAC  map  position  could  be  obtained  for  Cali  and  Buena  Ventura. 
The  position  at  Buena  Ventura  given  in  Table  9 is  the  one  Thompson  determined 
in  the  field  at  Buena  Ventura  from  navigation  chart  HO  1781  (Woollard  and 
Thompson's  Table  1,  p.  34).  It  does  not  seem  sufficiently  accurate  for 
comparison  of  the  map  offsets,  nor  consistent  with  the  differences  found  for 
the  other  sites. 

As  Table  9 shows,  there  exists  a systematic  offset  between  IGAC  map- 
and  satellite-determined  latitudes  and  longitudes;  so  the  satellite  fix  is 
between  0.100  and  0.180  min  to  the  south  and  between  0.150  and  0.190  min 
to  the  east  of  the  map  position.  The  offset  is  well  outside  the  error 
limits  that  can  be  assigned  either  to  the  map  or  to  the  satellite  position 
accuracy,  and  probably  relates  to  the  Bogota  datum. 

Height  relation  to  gravity 

Although  Woollard  and  Thompson  showed  that  anomalous  gravity  could  be 
related  qualitatively  to  uncertainties  in  position,  no  quantitative  relations 
were  established.  In  the  present  study,  and  excluding  the  data  for 
satellite  no.  64,  the  difference  between  map  height  and  satellite -determined 
height  (the  latter  height  including  geoid  undulations)  was  found  to  show  a 
systematic  relation  to  the  regional  simple  Bouguer  gravity  anomaly  of  about 
2 mgal  per  m height  difference  (see  Table  9,  last  two  columns,  and  Fig.  14), 
zero  Bouguer  anomaly  corresponding  to  w 32  m which  is  close  approximation 
of  the  geoidal  height  where  it  reaches  a maximum  value  in  western  Colombia. 


RELIABILITY  OF  POSITION  DATA  FOR  LONGITUDE  AND  HEIGHT 

SUVA  1971 


An  Important  question  to  be  answered  was  whether  or  not  the  Colombia 
longitude  and  heights  were  as  accurate  as  the  Internal  consistency  of  the 
data  suggested  by  their  cfy. 


Conventional  arithmetic  mean  longitude  and  latitude 

Since  many  longer  series  of  observations  were  reported  by  Daugherty, 
any  particular  series  at  a fixed  position  could  have  been  analyzed  by  the 
double  pass  method.  However,  C.  Marsh  (Hawaii  Institute  of  Geophysics) 
suggested  using  the  Suva  data  because  he  found  that  two  different  heights 
had  been  initialized.  This  circumstance,  not  reported  by  Daugherty,  promised 
to  provide  a particularly  good  opportunity  to  check  the  elevation  correction. 
In  addition,  a considerable  number  of  passes  had  been  recorded  that  were 
not  reported  by  Daugherty,  especially  during  the  first  few  days,  where 
initialized  antenna  height  (including  the  geoid  height)  was  75  m.  All 
available  passes  are  presented  in  Table  10A,  and  all  available  passes  where 
the  initialized  height  was  54  m are  listed  in  Table  11A.  According  to 
Daugherty,  the  information  contained  in  the  standard  fix  output  of  the 
HP2114  computer  used  was  based  on  the  MAPS-H-70356  program  aboard  the 
R/V  KANA  KEOKI  in  1971.  The  parameters  relating  to  the  movement  of  the 
receiver  were  deleted  and  the  format  of  the  remaining  quantities  was  changed 
as  shown  below: 

(1)  DATE  (DATE):  The  Julian  Day  Number  representing  the  consecutive 

numbering  of  the  days  in  the  year; 

(2)  TIME  (GMT):  The  Greenwich  Mean  Time  of  the  position  fix  in 

hours  and  minutes; 

(3)  SATELLITES  (SAT):  Five  near-polar-orbit  navigation  satellites 

were  available  for  observation.  A two-digit  numbering  code  for 
satellites,  related  to  the  orbital  semi-major  axis  in  kilometers, 
was  adopted.  The  correspondence  to  the  satellite  number  is 
indicated  below. 

SAT 

42 

54 

63 

64 

65 


SATELLITE  NUMBER 

30120 

30140 

30180 

30130 

30190 


12. 


The  numbering  code  (left  column)  is  higher  by  one  unit  than  that 
used  for  the  Colombia  data  by  Wool  lard  and  Thompson,  so  that 
SAT  42  (Suva  1971)  = SAT  41  (Colombia  1973)  and  so  forth. 

(4)  ELEVATION  ANGLE  (ELEV) : Vertical  angle,  measured  In  degrees,  of 

the  satellite  at  the  point  of  closest  approach  to  the  observing 
station; 

(5)  GEOMETRY  (GEOM) : Indicates  the  satellite  travel  is  either 

North  (N)  or  South  (S)  bound,  to  the  East  (E)  or  West  (W)  of 
the  observer; 

(6)  LATITUDE  (LATITUDE):  Latitude  of  the  observation  point  on  the 

satellite  datum  measured  in  degrees,  minutes,  and  seconds  from 
the  equator; 

(7)  LONGITUDE  (LONGITUDE):  Longitude  of  the  observation  point  on  the 

satellite  datum  measured  in  degrees,  minutes,  and  seconds  from 
the  Greenwich  Meridian  to  the  station; 

(8)  ITERATIONS  (IT):  Ntimber  of  iterations  required  for  the  program 

to  converge  on  the  positional  fix; 

(9)  DOPPLER  COUNTS  (CTS) : Number  of  24-sec  interval  counts  used  in 

the  computation  of  the  position  fix; 

(10)  DOPPLER  COUNT  SEQUENCE  (CT9Q):  Number  of  balanced  (symmetric) 

24-sec  counts  about  the  point  of  closest  approach. 

(11)  The  arithmetic  mean  position  in  longitude  and  latitude  was 
calculated  for  each  of  the  two  tables  separately  by  retaining 
only  those  passes  for  which  the  following  criteria  were  met: 


(1) 

15°  < ELEV  £ 75° 

elevation  angle 

(2) 

IT  < 4 

iterations 

(3) 

Dev  < 10. 0 arc  sec 

position  deviation 

The  deviations  from  the  mean  latitude  and  longitude  are  also  given 
in  the  Tables  10A,  11A,  as  well  as  the  rejection  criteria  used 
(last  three  columns).  The  mean  latitude  and  longitude  together 
with  their  statistical  reliability  are  given  in  Tables  10B 
and  11B,  and  the  mean  values  for  each  satellite  separately  were 
presented  in  Tables  10C  and  llC. 


13. 


Double  pass--longltude  and  height  results 

All  double  passes  (e.g.  two  consecutive  passes  of  the  same  satellite) 
in  Table  10  and  11  have  been  analyzed  for  longitudinal  position  and 
receiver  height.  Figures  8 through  13  show  the  f(K^)  versus  data  for 
12  consecutive  1971  days  (203  to  214)  for  Suva.  H 75  or  H 54  indicates 
the  initialized  antenna  height.  Satellites  are  identified  by  their 
number  (i.e.,  42),  direction  of  travel  (S  or  N) , and  the  chronological 
order  of  the  center  time  of  the  two  consecutive  passes:  1 or  2--to  9 

(on  day  204).  The  longitude  determined  from  the  double  pass  is  the  position 
where  the  connecting  line  between  the  two  passes  crosses  the  x-axis;  the 
inverse  slope  of  the  line,  corrected  for  latitude,  gives  the  antenna  height 
correction. 

It  is  clear  from  the  crossings  of  the  x-axis  that  one  can  easily 
distinguish  good  days  (like  203  and  209)  from  bad  days  (like  210),  and  that 
bad  double  pass  positions  like  63S5  (at  178°25".5800E)  on  day  203  can  be 
rejected.  Tables  12  and  13  give  longitudes  and  heights  (e.g.  H + correction) 
for  all  retained  double  passes,  averages,  standard  deviations,  and  standard 
deviations  of  the  mean,  ordered  by  N-  or  S-bound  satellites  or  by  days,  and 
for  all  satellites  and  days  taken  together. 

It  should  be  noted  that  for  Suva,  E is  in  the  positive  X direction  on 
Figures  11  through  16,  whereas  E is  in  the  negative  X direction  for 
Figures  1 through  7 pertaining  to  the  Colombia  longitudes.  As  a consequence, 
the  antenna  height  correction  changes  sign  between  the  two  sets  of  figures 
for  the  same  apparent  slope. 

Figures  15  through  18  are  graphical  summaries  of  the  contents  of 
Table  12  and  13.  The  deviation  of  the  daily  mean  longitude  from  the  total 
mean  value  only  once  exceeds  8 m (Table  12,  Fig.  15)  whereas  the  daily 
mean  heights  show  somewhat  larger  variations  (Table  13,  Fig.  16),  mainly 
due  to  the  lower  height  determined  from  the  mean  of  satellite  65.  This  is 
the  same  satellite  (no.  64  Colombia)  that  showed  a much  lower  height  for 
Laguna  la  Cocha  (Fig.  3)  two  years  later.  The  lower  height  is  also 
clearly  borne  out  in  Figure  16,  where  not  even  one  height  determined  from 
no.  65  (Suva)  (identified  by  +)  reaches  the  mean,  as  determined  by  all 
double  passes.  The  mean  height  determined  from  no.  65  of  32.3  ± 3.5  m 
deviated  by  more  than  26  m from  the  mean  of  the  other  four  satellites  of 
58.7  ± 2.4  m.  The  cause  of  this  rather  significant  discrepancy  is  not 
known,  but  as  Table  12  and  Figure  18  indicate,  it  does  not  seem  to  affect 
the  longitude  observations.  In  Figure  18,  the  mean  longitude  and  the  mean 
height  together  with  the  2a  error-length  are  indicated  for  each  satellite, 
and  the  positions  marked  S-  or  N-  are  the  corresponding  values  for  the 
southbound  or  the  northbound  orbits.  The  mean  of  all  the  satellite 
longitudes  and  the  mean  height  without  satellite  no.  65  are  indicated  by 
the  dashed  lines  together  with  the  2a  error  of  the  mean. 
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Finally,  the  writer  attempted  to  detect  any  possible  variations  as  a 
function  of  time  of  day  that  might  depend  on  changes  in  ionospheric  or 
tropospheric  conditions.  In  Figure  17  the  64  double  pass  longitudes  are 
shown  as  a function  of  the  time  of  day.  Values  that  have  been  rejected 
from  visual  inspection  on  a particular  date  are  surrounded  by  a square. 

The  two  rejected  double  passes,  closer  to  the  mean  than  some  of  the  non- 
re  jected  ones,  both  clearly  were  outside  the  remaining  seven  positions 
on  day  204  (Fig.  8,  Sat  64S2  and  64N7)  and  probably  would  have  been 
retained  on  another  date.  The  data  seem  to  show  a perhaps  slightly  tighter 
clustering  between  14  and  18  hours  GMT,  which  are  the  local  early  morning 
hours  (from  2 to  6 a.m.). 
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STANDARD  ERROR- -COMPARISON 


In  the  following,  discussion  will  be  limited  to  longitude  data,  since 
no  attempt  was  made  to  correct  the  latitude  data.  As  was  to  be  expected, 
the  accuracy  of  the  double  pass  longitude  determination  is  considerably 
higher  than  that  of  the  routinely  used  arithmetic  averaging  especially  in 
the  case  where  the  antenna  height  was  incorrectly  initialized.  The  extreme 
example  is  represented  by  the  longitudes  determined  for  Laguna  la  Cocha 
(Table  7),  where  the  a for  a single  observation  is  1'.41  (min  of  arc)  when 
the  arithmetic  mean  is  taken,  and  only  O'. 0047  for  the  double  pass  method. 
Even  the  more  realistic  comparison  of  the  ct  of  the  mean  yields  0.'408  for 
the  former  versus  O'. 0033  for  the  latter  method,  approximately  two  orders 
of  magnitude  better. 


Arithmetic  mean  versus  double  pass  method 

Returning  to  the  Suva  data,  we  see  that  the  set  with  the  eironeously 
initialized  height  of  75  m shows  a standard  deviation  for  a single 
observation  of  (1.9  sec  or)  O'. 032  as  compared  to  (1.5  sec  or)  O'. 025 
(Tables  10B,  llB)  for  the  set  with  the  nearly  correct  antenna  height  of 
54  m.  In  order  to  make  a comparison  with  the  double  pass  method,  the 
preceding  values  should  be  divided  by  /I  since  in  the  latter  method 
each  position  is  determined  from  two  consecutive  passes.  Therefore,  even 
in  the  best  possible  condition  of  nearly  correct  antenna  height  initialization, 
the  cr  (two  observations)  is  O'. 025 //2  .=  O'. 018  as  compared  to  O'. Oil  for  the 
double  pass  method  (Table  12).  The  o of  the  mean  also  would  be  comparably 
smaller  since  a total  of  137  passes  were  accepted  for  the  arithmetic  mean 
method  (see  "NSD"  in  Table  10B  and  in  Table  11B)  versus  64  double  passes 
(128  single  passes).  The  improvement  is  about  1.6.  However,  as  Table  10B 
shows  (for  antenna  height  error  of  about  20  m) , the  a (single  observation) 
increases  to  (1.9  sec  or)  O'. 032  or  0'.032//2  = O'. 022,  which  results 
in  an  improvement  of  2.0  for  the  double  pass  method. 

Double  pass  method:  fixed  shipboard  stations  versus  fixed  land  stations 

The  difference  in  accuracy  between  the  shipboard  data  and  the  land- 
based  data  (in  longitude)  is  somewhat  more  difficult  to  evaluate,  since  no 
single  common  set  of  passes  is  available  and  data  from  a single  land 
station  in  Colombia  have  at  the  most  three  double  passes.  As  discussed 
earlier,  the  results  of  all  accepted  double  passes  for  Suva  (Table  12) 
indicate  a a (single  ship  observation)  of  O'. 0111  or  19.6  m.  However,  in 
addition,  some  double  pass  positions  were  rejected.  To  evaluate  the 
standard  deviation  for  the  land  sites  the  following  expression  was  used: 
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where  AX  is  the  deviation  of  the  i**1  station's  j1"  observations  from 
the  mean  position  of  the  i*"*1  station;  and  n ■ 12  is  the  total  number  of 
double  pass  positions  obtained  (for  Buga,  Cali,  Popayan  I,  Laguna  la  Cocha 
and  Buena  Ventura;  see  Figs.  3 through  7).  The  value  thus  obtained  of 
ct  (single  land  observ.)  = O'. 0045  (or  8.35  m)  is  less  than  half  of  that  for 
the  shipboard  data.  It  should  be  emphasized  strongly  that  none  of  the 
double  passes  recorded  in  Colombia  had  to  be  rejected. 

Possible  explanations  of  the  difference  in  accuracy  are: 

(1)  ship  motion  near  the  docking  site;  and 

(2)  the  satellite  orbital  programming  and  the  navigator  programs 
used  in  1973  in  Colombia  were  more  accurate  than  those  used 
during  1971  in  Suva. 

I 
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CONCLUSIONS 


(1)  Except  for  maps  from  the  Buena  Ventura  area,  the  IGAC  and  TEXACO 
maps  agree  with  sufficient  accuracy  to  allow  the  determination  of  the 
seismic  refraction  stations  to  a common  datum  during  project  Narino  in  1973. 

(2)  For  all  Colombia  sites  for  which  IGAC  map  coordinates  are 
available,  the  satellite  navigation  positions  are  to  the  South  (by  O'. 103 
to  O'. 180)  and  East  (by  O'. 152  to  O'. 191)  of  the  map  positions.  The 
combined  (map  position  and  satellite-navigator  fix)  longitude  error  should 
be  smaller  than  30  m and  could  be  improved  to  less  than  10  m if  the  sites 
were  surveyed.  The  standard  longitude  deviation  for  a single  double  pass 
observation  was  8.4  m. 

Due  to  some  error,  probably  internal,  of  the  satellite  navigator 
computer,  elevations  are  reliable  only  for  Buga,  Tumaco,  and  Buena  Ventura 
and  are  probably  reliable  with  20  m or  less.  Elevations  for  the  remaining 
three  sites  should  be  used  with  caution. 

(3)  The  sliipborne,  fixed-point  satellite  navigation  determinations  of 
longitude  and  height  show  a standard  deviation  for  a single  double  pass 
observation  of  less  than  20  ra  for  all  five  operational  satellites  combined. 
However,  significant  differences  for  the  height  have  been  observed  for 
satellite  65  when  compared  with  the  others.  The  improvement  in  longitude 
position  is  at  least  a factor  of  1.6  over  that  obtained  by  the  usual 
arithmetic  averaging  and  increases  with  the  antenna  height  error.  For 
example,  for  the  Suva  longitude  the  improvement  was  a factor  of  2 for  an 
antenna  height  error  of  around  20  m.  The  standard  deviation  of  the  mean  in 
Suva  was  less  than  2.5  m in  either  longitude  or  height  over  a total  of 

12  days  of  recording  time. 

This  accuracy  suggests  that  many  more  harbor  sites  can  be  located 
with  respect  to  the  satellite  navigation  datum  within  8 m for  each  coordinate 
from  a single  day's  observations,  and  the  local  geoid  height  can  be  obtained 
as  well,  both  a considerable  advantage  over  the  arithmetic  averaging  method. 

(4)  Differences  in  accuracy  are  noticed  among  individual  satellites 

(at  least  for  the  Suva  data).  Whereas  the  o (single)  = O'. 0111  for  all 
satellites  combined,  the  a (single)  for  no.  64  was  only  O'. 0061  and  that  for 
no.  63  was  O'. 0092.  Also,  lower  deviations  were  observed  for  northbound 
than  for  southbound  orbits  of  individual  satellites  (exception:  no.  42). 

In  the  extreme  case  of  no.  64,  cr  (single  N-bound)  = O'. 0037  versus  CT  (single 
S-bound)  = O'. 0089.  One  observation  of  concern  was  that  satellite  65 
indicated  an  antenna  height  lower  by  more  than  26  m than  those  of  the  other 
satellites  used. 

In  conclusion  it  appears  that  the  satellite  data  can  be  used  to 
determine  geodetic  coordinates  (longitude,  latitude  and  height)  on  a 


' v' 


18. 


worldwide  basis  with  a reliability  equal  or  close  to  that  obtained  by  other 
more  costly  surveying  methods.  Bender  (1974),  for  example  indicates  that 
there  are  very  few  locations  in  the  world,  and  these  limited  to  continental 
areas,  where  geodetic  coordinates  have  been  established  with  a reliability 
of  2 m. 

Finally  it  should  be  pointed  out  that  the  question  of  the  corrections 
applied  and  this  report  itself  probably  never  would  have  been  made  if  the 
Magnavox  706  satellite  navigator  had  calculated  the  fixes  with  the 
initialized  antenna  heights  instead  of  what  it  probably  did:  used  a 

height  that  resulted  from  upward  shift  of  one  binary  digit  at  elevations 
above  1024  m. 
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Addendum 

After  this  report  was  completed,  Mr.  Wlllian  Unger,  University  of 
Wisconsin,  having  received  a preliminary  copy,  commented  on  (written 
communication,  October  24,  1975)  and  sent  calculations  transforming  the 
map  coordinates  from  the  Bogota  datum  to  those  for  the  NWL-9D*  coordinate 
system  used. 

The  following  table  gives  the  resulting  differences  in  latitude  and 
longitude  (in  minutes  of  arc),  as  calculated  by  Unger  from: 


Map  Coordinates  from  Table  9 for  Stations  Buga,  Popayan  I, 
Laguna  la  Cocha  and  Tumaco 


NWL-9D  Minus  SAT-NAV  Positions 
(in  minutes  of  arc) 


A LAT 

Buga 

Popayan  I 
Laguna  la  Cocha 
Tumaco 


+ 0.0071 

- 0.0273 

- 0.0720 
+ 0.0061 


A LONG 

- 0.0014 
+ 0.0077 

- 0.0298 

- 0.0089 


Besides  the  Laguna  la  Cocha  differences,  for  which  the  map  positions 
could  not  be  established  as  well  as  for  the  other  stations  (error  estimates 
with  ? in  Table  9),  only  the  Popayan  I latitude  is  outside  the  combined 
estimated  map  position  error  plus  the  Oj^  of  the  SAT-NAV  determination. 
However,  if  the  arithmetic  mean  value  of  the  SAT-NAV  latitude  is  used 
(Wool lard  and  Thompson's  Table  4,  p.  48)  the  A latitude  for  Popayan  I is 
- 0.0068  and  within  the  combined  errors. 


*The  NWL-9D  system  is  the  one  used  for  the  Transit  Navy  Navigation 
Satellite  System.  (See  Offshore  Technology  Conference  paper  OTC  1789 
by  Thomas  A.  Stansell,  1973.) 
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Table  1 


Site 

Passes 

Std.  Deviation, 

Tec 

Std.  Deviation  of  Mean, 

sec 

llonolu  1 u 
Lat 

79 

1.67 

0.19 

Long 

1.59 

0.18 

Pago  Pago 
Lat 

24 

1.56 

0.32 

Long 

1.63 

0.33 

Suva  (1) 
Lat 

^9 

1.86 

0.35 

Long 

2.53 

0.47 

Rabaul 

Lat 

6 

2.25 

0.92 

Long 

2.43 

0.99 

Guam 

Lat 

13 

1.88 

0.52 

Long 

1.32 

0.37 

Ponape 

Lat 

25 

1.93 

0.39 

Long 

1.62 

0.32 

Palau 

Lat 

48 

1.16 

0.17 

Long 

1.84 

0.26 

Suva  (2) 
La  t 

153 

1 .84 

0.15 

Long 

1.76 

0.14 

Wellington 

Lat 

46 

0.73 

0.11 

Long 

1.13 

0.17 

Table  2 


Sill.' 

L lev  . , 

m 

8 Pa  Bui- u 

W Fa sues 

K 

vs  W, 

sec 

Buna 

968 

3 

3 

Lai 

Long 

1.014N 

2.418E 

Call 

1326 

3 

3 

Lat 

Long 

0.36S 

63.78W 

Popayan  1 

1743 

2 

2 

Lat 

Long 

0.66N 
48  - 51W 

Popayan  It 

7 

1 

2 

Lat 

Long 

3. 60S 
61.98W 

Laguna  da  Cocha 

3133  * 

5 

5 

Lat 

Long 

2.73S 
155. 27W 

Tumaco 

5 

3 

3 

Lat 

Long 

3.902N 

0.36W 

Buena  Ventura 

5* 

3 

3 

Lat 

Long 

2 . 142S 
0.462W 

Bogota 

2650 

1 

4 

Lat 

Long 

7.72N 
278. 01W 

Table  taken  from  Wool  lard  and  Thompson,  p.  38. 

* These  elevations  are  In  error  In  Wool  lard  and  Thompson. 
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TABLE  3 


ILXACU  MAP  X,  Y.  COORDINATES  CONVERTED  TO  GEOGRAPHICAL  COORDINATES 
LSY  COMPUTER  PROGRAM 


CM  1 I -III  IN  79-1  5W 

x 6 O 0.0  000 

f :L  - 1 ) -005  78-  TOW 

602  3 10.0000 
f : I - 2 l-OON  7/-V5W 

1 1 0 r'.  r 20.0000 

f : l - 3 l-OON  77-OOW 
(.02  1 60.0000 
f:t-9  l-OON  76- 1 6W 
6021 10.0000 
C : A-OON  77-OOW 

Dpi  320.0000 
r:J-'.  3 00N  77-OOW 
ft?  i<530. 0000 


Y 5?  i860. 00 00 


N 0°  So'.N/OS  VV  /9<"  l<i'.V989 


60  7 7 10.0000 

66 


0 69. 9769 


>'  I OiO. 00OO 


I)  59.9831 


( 7't'JSO.OOOO 


0 99.5960 


7 6 86  NO. 0 000 


1 0.0061 


75690. OOOO 


3 59.9930 


6 75320.0000 


3 0.0110 


7 8 2V.W-»‘-o 

. 992 


7 7 66.991'’ 


76  69.9968 


76  15.0096 


76  59.9961 


77  0.0058 


TABLE  4 


COMPARISON  OF  REFERENCE  POINT  POSITIONS  ON  IGAC  AND  TEXACO  MAPS 


SUC^IA  Vt«T'JftA-VAI  1 c 

260-1 ll-t  J-4 

979520.0000 

1011400.0000 

3 

99.9307 

76 

51.6424 

924S0O.OOOO 

670060.0000 

3 

94.9363 

76 

90.6497 

0 

0.9944 

0 

-0.0573 

260-1 Il-C  J-4 

424700.0000 

1012330.0000 

3 

99.1292 

76 

60.1941 

929200.0000 

670930.0000 

3 

94.9193 

76 

90.2209 

0 

0.6139 

0 

-0.0340 

260- I ll-C  J-4 

926230.0000 

1013310.0000 

3 

99.9190 

76 

97.6646 

929620.0000 

679910.0000 

3 

99.2866 

76 

97.7020 

0 

0.6293 

0 

-0.0302 

260-lll-C  J-4 

920630.0000 

1012470.0000 

j 

92.077* 

76 

40.1103 

921020.0000 

677620.0000 

3 

92.70S6 

76 

90.9290 

0 

0.0070 

0 

(1.0110 

260- I Il-C  J-4 

921920.0000 

1013390.0000 

3 

53.3602 

76 

97.621  7 

921  100.0000 

6 70600.0000 

3 

53. 2034 

76 

58.3593 

0 

0. 1568 

0 

-0.7376 

279-1 -A  J-4 

919760.0000 

1013690.0000 

3 

92.4161 

76 

57.4594 

920180.0000 

679120.0000 

a 

92.3372 

76 

68.1189 

0 

0.0709 

0 

-0.6592 

MIHIRA-VAILF 

20O-II1-D  J-4 

089940.0000 

1009970.0000 

3 

39.9072 

76 

16.4069 

009760.0000 

796  180.0000 

3 

39.  9604 

76 

16.4421 

0 

0.0102 

0 

-0.0096 

290-111-0  J-4 

007060.0^00 

1O0S4AO.OOOO 

3 

39.0776 

76 

10.6957 

000000.0000 

792140.0000 

3 

39.0520 

7i» 

10.6699 

0 

0.0249 

0 

0.0259 

20O-III-O  J-4 

•03970.0000 

1066190.0000 

3 

32.9660 

76 

18.3141 

0041  90.  (1000 

752  799.0000 

3 

32.9445 

76 

1H.J195 

0 

0.0223 

0 

-0.009<r 

200-1 fl-0  J-4 

000970.0000 

1O05H3O.OOOO 

3 

31.1223 

76 

10.5100 

000/60.0000 

792400.0000 

3 

31.0044 

76 

19.4770 

o 

0.0379 

O 

0.0331 

(|-OI»|  . ) 

TABLE  4 (coni.) 


Mf  ACAOFPf 

096600.0000 
44S6lO.O<HjO 

6 66O00.l)0i»O 
68o?0O.'ii>ftt> 

6 9 | 100. 00O0 
oSi 1020. OOrtO 

69rt«»5rt.O0t>0 

6909io.*»noc 

'•91?  90.00*0 
69)560.0000 

TUQUf*f!-N8»lN0 

60^20.0000 

610150.0000 

609520.0000 

605790.0000 

60)600.0000 

60)700.0000 

608110.0000 
608)50.0000 

603)10.0000 

60)620.0000 

60)660.000(1 

606180.0000 


)Ht»-J-l*  1-1 

6*11*100.0090 
6576(0.0000 

386-IV-O  L -3 

691 700.0000 
667  8 7(1.0000 


)8  6- 1 V-8  t-3 

<>602  70.0000 

656660.0000 

336-lV-S  1-) 

990660. 00*10 

656690.0000 

)86- 1 V-8  1-) 

999610.0000 
65561  '1.O0U0 

629-1 V-C  l-f 

923190.0000 

589090.0000 

626- I V-C  1-2 
922720*0000 

598600.0000 

628-IV-C  L *2 

926000.0000 

589860.0000 

628-IV-C  1-2 

926310.0000 

690190.0000 

62  6-1 V— 0 1-) 

931660.0000 
5975)0.0000 

628-tV-n  l-l 

936830.0000 

600670.0000 


1 65. )560 

l 65.3116 

0 0.0)96 

1 65.5671 
l 66.5720 

0 0.0663 

1 68.3363 

1 68.13)0 

0 18.2811 

1 68.0902 
l 68.0737 

0 0.0165 

1 65.5172 
l 69.5055 

0 0.0087 

1 9.2905 

1 6.2756 

0 0.0069 

1 6.06)6 

1 6.0750 

0 -0.0116 

l 0.76)0 
l 0.7835 

0 -0.0605 

1 3.2986 

1 9.2909 

0 0.0077 

1 0.6950 
I 0. 7669 

0 -0.0699 

1 1.0)71 

1 1.0501 

0 -0.0129 


77  9.222? 

77  9.1659 

0 0.0273 

77  9.2879 

77  9.2696 

U 0.0379 

77  10.1018 

77  10.0239 

0 V0.0779* 

77  9.5915 

77  9.8892 

0 0.002) 

77  10.6577 

77  10.6731 

0 -0.0156 

77  66.Z594 

77  66.2122 

0 0.0672 

77  66.5127 

77  66.6702 

0 0.0625 

77  65.8220 

77  65.799) 

0 0.0227 

77  65.6556 

77  65.61)9 

0 0.0615 

77  61.6769 

77  61.6624 

0 0.0145 

77  39.9645 

77  39.97)9 

o o.ono 

(COflt.) 
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TABLE  4 (coni.) 


t 


629-III-6  1-3 


M | 994.0000 

*♦5*1360.0600 

1 

5.60  71 

77 

81721.1.0040 

(•ZlMiZO.oflno 

1 

5*6730 

n 

H.ttl) 

0 

-0*01-6 

0 

-0.017. 

628-111-0  L-J 

9&707iI*0000 

1 

8* |46h 

77 

*2.121* 

61  y\?  11.0000 

63  3900* 0000 

l 

8.  Id  76 

77 

22.1003 

0 

0.0>9? 

0 

0.0707 

ClUSTQ-MiRlNO 

629-1 i-r  i - i 

825250.00OO 

900610.0000 

l 

12*6031 

77 

15.7001 

*25470.0000 

64G690.O000 

1 

12*  6106 

77 

IS. 2221 

-0.0075 

0 

-0.0220 

629— 1 1— C 1-3 

624950.0000 

98J400.0000 

1 

12*6606 

77 

13.(039 

6251  10.0000 

669320.0000 

1 

12.6173 

77 

13.6062 

u 

0.0272 

0 

-0.0926 

429-ll-C  1-3 

628130*0000 

° 77190*0000 

1 

14*1657 

77 

17.1*17 

628360*0000 

643090.  0000 

1 

16*2740 

77 

17.1670 

0 

-0.0083 

0 

-0.0153 

425-11-C  t-3 

624690*0000 

982670.0090 

L 

12*  1420 

77 

16.3052 

62457  0*0003 

668370*0000 

l 

12.1241 

77 

16.3173 

0 

0.  31 79 

0 

-0.0121 

CJNfCOA-PUTunAYO* 

630-1 V-C  1-6 

602600.0000 

m 1060630.0000 

1 

0.3120 

76 

67.0616 

602610*0000 

▼06320.0000 

l 

0.2566 

▼ 6 

63.1071 

0 

0.0577 

0 

-0.065? 

630-  IV- 0 1.-6 

603?  20*0000 

1068000.0000 

1 

0.7076 

76 

34.9806 

603340*0000 

7139.70  . 700  0 

1 

0*6533 

76 

39.0256 

0 

0.0543 

0 

-0.0666 

430—1 V—D  1-6 

607970.0000 

1049260*0000 

1 

3.725? 

76 

30.3013 

607570*0000 

715090. 700 0 

1 

3*1636 

76 

34.3664 

o 

0.0616 

0 

-0.0453 

430-TV-D  1-6 

6 06850. 4900 

1047630.0000 

1 

2*6176 

76 

34.1400 

606  800*0000 

711500.0000 

1 

2.5287 

76 

39.7673 

0 

0.0009 

0 

-0.(3625 

ClPTO-lSIS-6uTUf14rO 

46  7-I-C  tt-6 

566980*00(10 

1063680*0000 

0 

31.2147 

76 

30.5635 

345000*0000 

729800.0000 

0 

31*2009 

74 

30.6663 

17 

0*0158 

0 

0.0972 

467-1 -C  W-4 

466610*0000 

1043970.0000 

0 

30.OJ02 

76 

30.3764 

564 430.0000 

730110*9000 

0 

10.0246 

76 

10.2790 

0 

0.0147 

0 

0.0476 

- - . . _ — 


Ti 


f 


TABLE  5 


Summary  of  Results  of  Table  4 by  Area 


IGAC  North 
of  Texaco 

IGAC  West 
of  Texaco 

Number  of 
Points 

Area 

Minutes  of 
Arc 

a 

Minutes  of 
Arc 

a 

Buenaventura-Valle 
North  of  922 

+ 

0.613 

0.015 

- 0.043 

0.010 

3 

South  of  922 

+ 

0. 108 

0.035 

- 0.736 

0.062 

3 

Palmira -Valle 

+ 

0.026 

0.007 

+ 0.020 

0.018 

4 

Mercaderes-Cauca 

+ 

0.027 

0.015 

+ 0.013 

0.021 

4 

Tuqueres-Narino 

- 

0.017 

0.021 

+ 0.030 

0.014 

6 

Yacanquer-Narino 

- 

0.004 

0.014 

+ 0.002 

0.019 

2 

Pasto-Narino 
(east  of) 

+ 

0.006 

0.014 

- 0.013 

0.007 

4 

Mocoa-Putamayo 
(west  of) 

+ 

0.066 

0.014 

- 0.060 

0.017 

4 

PTO  As i a -Put amayo 

+ 

0.015 

0.001 

+ 0.097 

0.000 

2 

Average  for  32  points 

0.018 

0.016 

32 

Average  for  29  points 

excluding  Buenaventura- 
South  of  922 

0.015 

0.015 

29 

(No  good  reference 
section) 

points 

were  available  In 

the  Buenaventura-South 

Coordinate  corners 

(see 

text) 

From  Texaco  map 

+ 

0.0082 

0.0130 

+ 0.0008 

0.0076 

7 

Table  6 


Error  In  Longitude  Position  Away  from  True 
Position  per  Meter  Error  In  Antenna  Height  with 
Increasing  Satellite  Elevation  Angle  , 

Based  on  Stansell  (1970) 

Error  f(K):  in  nautical  miles  per  meter  error  in 

antenna  height 

Elevation  K:  degrees 


K (degrees) 

f(K)  10“^  n.  miles 

10 

2.81 

20 

3.62 

30 

4.62 

40 

6.05 

50 

8.24 

60 

11.63 

70 

17.55 

'•  i f ~ nil  ' ■ 7 


A 

J 


■ 

Satellite 

Navlgat Ion 
(from 

TABLE  7 

Results  for  Each  Observation 
Woollard  and  Thompson) 

Site 

Initialized 

Height* 

Elev. 

» 

Ceold 

Antenna 

Time 

Sat. 

Dir 

K 

Lat . 

Long . 

8 

H 

Buga,  Colombia; 

Feb. 

12,  1973; 

elev 

. 968  m 

8 

960 

0700 

41 

S 

46E 

3*53.654' 

76*18.137' 

0850 

41 

S 

16W 

53.682' 

18.200' 

0946 

53 

S 

38E 

53.655' 

18.158' 

1034 

64 

N 

24E 

53.621' 

18.175' 

1136 

53 

S 

18W 

53.679' 

18.198' 

1218 

64 

N 

37W 

53.648' 

18.193' 

1418 

62 

N 

64W 

53.632' 

18.197' 

Call,  Colombia; 

Feb. 

14,  1973; 

elev 

. 1326  m 

8 

1326 

0252 

62 

S 

2 9W 

3*25.725' 

76*34.018' 

0538 

63 

S 

16E 

25.721' 

34.768' 

0700 

41 

S 

49E 

25.783' 

35.224' 

0726 

63 

S 

44W 

25.700’ 

33.798' 

1042 

64 

N 

35E 

25.673' 

35.027' 

1228 

64 

N 

26W 

25.770' 

34.014’ 

Popayan 

(1),  Colombia;  Feb.  14 

'15, 

1973;  elev. 

1743  m 

8 

1743 

2130 

53 

N 

29E 

2*26.585' 

76*37.043' 

2215 

64 

S 

15E 

26.695’ 

36.863' 

2316 

53 

N 

2 5W 

26.647’ 

36 .115' 

0004 

64 

S 

3 7W 

26.611' 

35.974' 

Laguna 

de  la  Cocha, 

Co lomb 1 a ; 

Feb . 

18,  19,  and 

20,  1973; 

15 

2808 

elev . 

3133m 

1918 

63 

N 

33W 

1*08.172' 

77*07.748' 

2012 

41 

N 

22W 

08.189' 

08.025' 

2134 

53 

N 

48E 

08.004’ 

10.440' 

2238 

64 

S 

37E 

08.242' 

09.897' 

i 

2322 

53 

N 

14W 

08.145* 

08.143' 

i. 

0026 

64 

S 

16H 

08.145' 

08.170' 

0216 

62 

S 

40W 

08.047' 

07.515* 

1822 

41 

S 

35E 

08.050' 

09.905' 

1928 

63 

S 

2 2W 

08.197  ' 

08.003  ' 

^ . 

2008 

41 

S 

20W 

08.175' 

08.058' 

* f 

2138 

53 

s 

6 4E 

07.92  ' 

11.659’ 

2248 

64 

N 

56  E 

08.312’ 

10.797' 

(cont.) 

Ilk 

Imrnamm 

. : „ 

1 1 > Ifct  1 — - 
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TABLE  7 (cont.) 


Time 

Sat  . 

Dir. 

E lev  . , 
K 

Lat . 

Long  . 

Initialized 

Height* 

Geold  Antenna 

8 H 

Tumaco , 

Col omb 1 a ; 

Feb. 

22,  1973; 

e lev . 5 m 

15 

5 

1818 

41 

N 

33E 

1*49.139' 

78*43.846' 

1940 

63 

N 

17W 

49.084', 

43.840' 

2142 

3 i 

N 

74E 

49.136' 

43.854' 

2258 

64 

S 

74F. 

49.17/ 

43.838' 

Buenaventura,  Colombia 

; Feb.  27, 

1973;  elev.  5 

m 

15 

10 

2100 

53 

N 

52E 

3*53.522' 

77*03.939' 

2226 

64 

S 

6 9E 

53.535’ 

03.903' 

2248 

53 

N 

13W 

53.610' 

03.899' 

0014 

64 

S 

7V 

53.545' 

03.904' 

0200 

62 

S 

26W 

53.528' 

03.903' 

0454 

63 

S 

2 3 E 

53.519' 

03.887 ' 

Bogota , 

Co lombi a ; 

Mar. 

2-4,  1973; 

Elev.  2650  m 

10 

2650 

1258 

62 

N 

56W 

4*37.847' 

76*0.969' 

0512 

63 

S 

66E 

37.732’ 

6.75’ 

0628 

41 

S 

5 9W 

37.401’ 

1.530’ 

0700 

63 

S 

9W 

37.537' 

3.214’ 

1358 

62 

N 

12W 

37.658' 

2.753’ 

Remarks : 

; The  lat. 

and 

long,  columns  have  been 

taken  directly 

from 

Thompson's  original  Satellite  Navieatlon  Log. 
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TABLE  8 


COLOMBIA  SITE  POSITION  IGAC  MAP 


BOGA  HOTEL  GUAOELAhARA  261-11 I-D-4 
ELtV:  973  METERS  X Y 

922430.0000  1086062.0000 

PCPAYANI  HOTEL  MUNASTfcWIC  1 :2000  PGP  AY  AN  5 
hit v:  1 73 J .5  METERS 

762C04.0000  t 052 12  0.0000 

laguna  lacOcha  HOTEL  429-IV-B 
ELEV:  2 BO  8 METERS 

(>1  7230.0000 


N 


3 53.  8331 


2 26.8020 


1 e.2516 


992390. 0000 

TUMACU  4 M NE  CE  BENCHMARK  103TS-1.1951  383-lV-A 

ELEV!  3 METERS  FROM  SEA  LEVEL 

692990.0000  816090.0000  l 49.3142 


w 

76  1 e. 3650 


76  36.7369 


77  8.9574 


78  4 4. C 126 


TABLE  4 


K.Ai  Map  Coordinates  ami  Satellite-Determined  Coordinates  I nr  Colombia 

LATITUDE  NORTH  LONGITUDE  WEST 


Error 

Man  Estimate, 

m 

SAT  NAV 
Wool  lard, 
Thompson 

A 

LAT 

Error 

Map  Estimate, 
hi 

SAT  NAV 
Double  Pass 
Method 

o 

M 

m 

A 

LONG 

Buga 

3°  53.8331 

+10 

53.653 

0.  18» 

76° 

18.366 

+20 

18.1820 

+1.8 

0. 18« 

(’alt 

3’ 

- - 

25. 709 

- 

- 

34.4583 

_3I0.0 

— 

Popayan  I 

2°  26.802ft 

i 2 

26.656 

0.147 

76° 

36. 7389 

+ 3 

36.5477 

+21.1 

o . ; *> ? 2 

Laguna  1 a Cocha 

l”  08.2516 

+30  ? 

08. 149 

0.  103 

77° 

08.9574 

+20  7 

08.8050 

+ 6.2 

0.1524 

Tumaco 

1°  49.3142 

+1  5 

49.1  14* 

0.180 

78 ** 

44.0126 

+1  5 

43.8429* 

+ 2.5 

0.1697 

Buenaventura 

3*  53.60* 

+ ? 

53 . 54 1 

-0.06 

77° 

03.72* 

1 ? 

03.9065 

+ 7.5 

-0.1865 

*Civen  by 
Wool lard  and 
Thompson  from 
Chart  HO  1781 

* arith- 
metic 
average 

*Given  by 
Woollard  and 
Thompson  from 
Chart  HO  1781 

*From  best 
line  of  no 
41,53,63 
8 i ngle 
patterns 

HEIGHT 


Ma  p , 
m 

Error 
Es t i ma  t e , 
m 

SAT  NAV 
Double  Pass 
H + G + Corr 
m oM 

Height 
Di f f erence, 
m 

SAT  ; 
W i thout 
H + G + 
m 

NAV 

Sat  no.  64 
Corr 

a 

M 

Height 
Di f f erence , 
m 

Approx imate 
Bouguer 
Anoma 1 ies. 
mgl 

Buga 

9 73 

+10 

1010 

+16 

37 

1021 

+i0 

48 

+ 25 

Call 

-- 

*1341.5 

+79 

1397 

Popayan  I 

1 731 

+ 3 

*1  726 

+45 

-5 

1694 

-37 

-135 

Laguna  la  Cocha 

2808 

+ 20 

*2819 

H+30 

11 

2771 

+24 

-37 

-175 

Tumaco 

3 

+ 2 

■f  1 9 . 8 

+ 2 

16.8 

"t  19.8 

+ 2 

16.8 

0 

Buenaventura 

7* 

+ 1 

9 

+26 

2 

-9 

-16 

- 75 

*Measured  *Assuming 

from  slide;  digital 

above  water  jump  ( see 

text  ) 

"J’From  single 
passes  without 
no.  64 


TABLE  10  A 


R/V  *ANA  K E0K1  1971  POSITIONAL  OATA,  SUVA,  Fiji 
«00«E0  TO  OOIPHINS  AT  N/E  KING’S  WHARF.  ANTE NNA 


DATE 

GMT 

SAT 

EISV 

GEOM 

i at  i roue 

LONGITUDE 

;o2 

3a  0 

6 1 

.?  6 

. - ► 

lb 

7 

46.42'- 

l/H 

25 

3 i.  24r- 

;d3 

31 2 

A ' 

7 7 

1 H 

7 

a-.,  it: 

178 

25 

26. 8?  £ 

20  3 

A 6 

6-» 

A 7 

S-. 

lb 

7 

A 6 . 001 

1 78 

25 

31. 6dr 

<03 

512 

aC 

8 

S-rt 

lb 

7 

A } . 6 A ‘ 

178 

25 

2 6 • 2 n r 

203 

6*1  A 

6-* 

.1 

S-w 

1 8 

7 

A6.  OCS 

! 78 

25 

32.68* 

203 

75  2 

oS 

7 it 

N-t 

lb 

7 

A 7 « 6 A c. 

1 7o 

? 5 

? 4.9ft" 

203 

62  0 

5a 

7 

S-‘ 

lb 

7 

6H.50S 

170 

25 

26.  221" 

2D1 

9 u 

o2 

1 3 

N-.‘ 

2 6 

7 

Ail.  30$ 

! 78 

25 

32.  0d» 

503 

9A*> 

66 

7 

N-W 

In 

7 

a2 • 8A ' 

170 

25 

29.  AOr- 

50  3 

100  2 

5A 

77 

S-W 

18 

7 

•♦6  • 9PS 

178 

25 

A5.66r 

21'? 

! .04 

6 i 

53 

N-W 

18 

7 

AM. 1 6S 

170 

25 

32. 887 

50  3 

l a3  a 

A 2 

35 

N-t 

lb 

7 

A 7. 5 2* 

1 78 

25 

31.0?* 

20  3 

15A2 

6 A 

18 

N-r 

1 b 

7 

A 8 • i 8 c 

1 70 

25 

3?. OA  r 

50  3 

1622 

42 

23 

N-W 

10 

7 

A 7 . 1 0 j 

178 

25 

32.  76F 

203 

1726 

6 A 

A5 

N-W 

18 

7 

47.58S 

178 

25 

33. 2AT 

503 

1 9 T o 

65 

65 

S— * 

1 8 

7 

Ao  . 3 8G 

1 78 

25 

2 9.20E 

503 

2108 

63 

10 

s-e 

18 

7 

A5.A2S 

178 

25 

3a. d6" 

5‘»3 

212a 

5*. 

52 

N-.- 

18 

7 

A3.66S 

17P 

25 

33. OOC 

So  3 

.252 

6? 

69 

S-w 

18 

7 

AA.«2C 

178 

25 

34.  a*.  : 

50  3 

2312 

5A 

17 

N-W 

18 

7 

A7.0AS 

178 

25 

32.765 

70A 

220 

A> 

25 

S-E 

18 

7 

A5. 1 8S 

178 

25 

3 3.  OOF 

50  A 

330 

6A 

13 

$-1 

10 

7 

49 .86* 

178 

25 

33.66  = 

50  A 

AU6 

A2 

32 

i- * 

10 

7 

A4.52S 

1 78 

? 5 

3?. Aft  7 

20  A 

51*. 

6 A 

57 

S-w 

1H 

7 

46. 7e* 

176 

25 

33.  36  r 

5 *'A 

7»A 

65 

2B 

N— 2 

18 

7 

56.90* 

170 

25 

30.90  c 

50A 

050 

65 

25 

N-W 

18 

7 

AO. i PS 

J 7 C 

2 5 

32. 80C 

204 

512 

5A 

35 

S-E 

lb 

7 

AA. 6 AS 

176 

25 

32. 5«F 

50  A 

1016 

63 

57 

N-w 

18 

7 

AP.A0”> 

176 

25 

30.7  8-J 

;oa 

: 1 Ou 

5A 

?3 

S- w 

\ 0 

7 

A6. 091 

17b 

25 

31.  3R" 

50a 

12'.'A 

63 

l 3 

N-w 

18 

7 

A5.00S 

178 

25 

32. loe 

20  a 

1 3A2 

A 2 

9 

N-C 

ie 

7 

52.96* 

) 78 

25 

M.56E 

20  A 

1526 

A? 

72 

N-W 

18 

7 

50.003 

170 

25 

37.  26E 

50A 

163  8 

fc  A 

72 

N-F 

18 

7 

A 7. AOS 

170 

25 

28. 625 

50A 

i A.  6 

6 A 

10 

N-W 

1 b 

7 

A 7.2  0* 

1 70 

25 

3I.38C 

50  A 

1 850 

65 

26 

S-t 

18 

7 

AS. 12S 

178 

25 

30. 70T 

2 ’JA 

2036 

5A 

17 

N— *. 

18 

7 

A 8 . l 0S 

178 

25 

31.02- 

50  A 

250A 

6> 

55 

1 8 

7 

A5.O0S 

178 

25 

31.02F 

50A 

22?  A 

5A 

58 

N-W 

10 

7 

AH. 90S 

170 

25 

33.A2r 

50  A 

2350 

63 

17 

s-h 

18 

7 

AA.82S 

176 

25 

32. 2? e 

205 

120 

42 

7 

s-t 

18 

7 

53.35$ 

178 

25 

3A. 56r 

50  5 

310 

A? 

76 

S-. 

10 

7 

AA. 40$ 

178 

25 

22.56" 

505 

A'  A 

6 A 

55 

s-r 

18 

7 

55. 90S 

l 78 

25 

31.029 

50  5 

50  0 

52 

7 

S-w 

18 

7 

51. 18S 

170 

25 

31.02C 

50  5 

61  A 

6A 

15 

S-w 

1 H 

7 

AA.  95  3 

l 78 

25 

?!  .ABC 

50  5 

b >2 

65 

69 

N-W 

1(1 

7 

A6.92S 

178 

? 5 

3 7.  7AC 

205 

622 

5A 

il 

S-f 

18 

7 

57.  2 OS 

1 7d 

25 

34.03C 

205 

926 

63 

21 

N-L 

lb 

7 

AH. 48$ 

170 

25 

31. 

505 

1006 

5A 

65 

S-W 

10 

7 

56.96: 

17P 

26 

33.36" 

505 

1 1 1 A 

6? 

37 

N-W 

: 8 

7 

AM. OOS 

1 78 

25 

? 3 • y 6 F 

505 

1 A3  2 

A2 

37 

N-r 

18 

7 

58.255 

178 

25 

30. AM ’ 

505 

1550 

6A 

28 

N--. 

18 

7 

A r , 0 6 '• 

178 

26 

31.56' 

505 

It  i6 

55 

21 

N-W 

18 

7 

A 7.  * A 

176 

25 

33.  OO*’ 

205 

l 73  o 

65 

30 

N-W 

18 

7 

4"  . 34' 

1?0 

25 

33. 2 A F 

205 

lflOA 

65 

0 

S-r 

10 

7 

A 6 . 78S 

* 70 

25 

3 ? . * . r 

505 

19A0 

65 

71 

S-w 

18 

7 

65. 77S 

1 78 

27 

29.  7ft  F 

205 

2116 

63 

16 

s-r 

18 

7 

Aft. O0S 

170 

25 

33.  2AC 

<u5 

21  3 A 

55 

65 

N-r 

l 8 

7 

A 7. 02  > 

! 78 

25 

?0.2AC 

205 

2 300 

63 

50 

s-w 

10 

7 

Aft.  2 ft  5 

178 

25 

33.  3(>r 

505 

2 320 

5A 

12 

N-W 

10 

7 

AO.08S 

17P 

25 

31.267 

506 

216 

52 

27 

S-E 

18 

7 

A A.  A 6*. 

178 

25 

32. 3AF 

506 

33H 

65 

?l 

s-c 

1H 

7 

A6. 56$ 

178 

25 

32.22C 

206 

A05 

5? 

29 

S-w 

18 

7 

45*46? 

! 78 

25 

? 3.241- 

506 

5 ? A 

6A 

39 

S-w 

18 

7 

A6.  62r 

178 

25 

3 A.  o?r- 

506 

71  A 

65 

52 

N-r 

1ft 

7 

Aft. 6?S 

1 70 

25 

30.1 ?r 

506 

902 

65 

16 

N-W 

10 

7 

AO. S8S 

178 

25 

31.  9?r 

206 

51  8 

5a 

5 5 

S-e 

18 

7 
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Fig.  8.  Left:  Double  pass  correction  for  the  Suva  harbor  position  on- 

board R/V  KANA  KEOKI  for  day  203  in  1971;  the  number  below  the 
satellite  identification  gives  the  chronological  order  in  which 
double  passes  occurred.  Note  that  63S5  clearly  is  outside  of  the 
remaining  five  double  pass  longitude  positions  and  therefore  has 
been  rejected.  The  spread  of  the  data  is  only  25  meters.  The 
standard  deviation  of  the  mean  of  the  five  double  pass  observa- 
tions is  0.0031  minutes  of  arc,  or  5.5  meters. 

Right:  Double  pass  positions  for  day  204. 


Fig.  9.  Left:  Satellite  elevation  angle  function  (right  side  of  line) 

as  function  of  the  satellite  elevation  angle  (left  side  of  line),  as 
taken  from  Stansell  (1970);  see  text. 

Right:  Double  pass  positions,  Suva  harbor,  day  205,  1971. 
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10.  Double  pass  positions  for  Suva  harbor  onboard  R/V  KANA  KEOKI 
or  days  206  and  207,  1971. 
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Fig.  12.  Double  pass  positions  for  Suva  harbor  onboard  R/V  KANA  KE0K1 
for  days  211  and  212,  1971. 
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Fig.  13.  Double  pass  positions  for  Suva  harbor  onboard  R/V  KANA  KEOKI 
for  days  213  and  214,  1971. 
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SUVA  ANTENNA  HEIGHT  DOUBLE  PASS  METHOD 


Fig.  17.  Suva  longitude  positions  as  a function  of  the  time  of  day 
(in  GMT).  Note  the  tighter  clustering  in  the  early  morning  (local) 
hours  between  14  and  1/  hours  GMT. 
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FIG. 17 


Fig.  18.  Suva  harbor  longitudes  and  heights  for  individual  satellites. 
I.rror  bars  indicate  two  standard  deviations  of  the  mean.  Positions 
marked  N-  or  S-  are  the  means  of  the  northbound  or  southbound  double 
pass  observations.  Note  the  considerably  lower  height  determined 
from  satellite  no.  65,  well  outside  observational  errors. 
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period  of  a week.  It  is  also  shown  that  the  method  allows  identification 
of  satellites  not  giving  reliable  data  and  permits  conclusions  to  be  reached 
as  to  the  cause  of  an  apparent  error  dependance  on  surface  elevation  noted 
above  1000  ra  in  the  Andean  region  of  Colombia,  South  America,  reported  by 
Woollard  and  Thompson  (1974).  An  additional  advantage  of  the  double  pass 
method  brought  out  is  that  it  allows  geoidal  heights  to  be  determined  for 
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